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ABSTRACT 


ARCH  models  are  used  widely  in  analyzing  economic  and  financial  tme  series  data.  Many 
tests  are  available  to  detect  the  presence  of  ARCH;  however,  there  is  no  acceptable  proce- 
dure available  for  testing  an  estimated  ARCH  model.  In  this  paper  we  develop  a  test  for 
a  linear  regression  model  with  ARCH  disturbances  using  the  framework  of  the  informa- 
tion matrix  (IM)  test.  For  the  ARCH  specification,  the  covariance  matrix  of  the  indicator 
vector  is  not  block  diagonal,  and  the  IM  test  is  turned  out  to  be  a  test  for  variation  in 
the  fourth  moment,  i.e,  a  test  for  heterokurtosis.  An  illustrative  example  is  provided  to 
demonstrate  the  usefulness  of  the  proposed  test. 

Key  Words:  Autoregressive  Conditional  Heteroskedasticity;  Information  Matrix  Test;  Dou- 
ble Length  Regression. 
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l.INTRODUCTION 


In  a  seminal  paper,  Engle  (19S2)  introduced  the  autoregressive  conditional  hetero- 
scedastic  (ARCH)  models.  These  models  are  now  very  popular  in  analyzing  financial  and 
economic  time  series  data  [for  a  recent  review,  see  Bera  and  Higgins  (1993)].  There  are 
many  procedures  available  to  detect  the  presence  of  ARCH.  However,  estimated  ARCH 
models  are  not,  in  general,  tested  thoroughly,  possibly  because  there  is  no  acceptable 
procedure  for  doing  that.  In  this  paper  we  derive  a  simple  specification  test  for  an  estimeted 
ARCH  model  in  the  linear  regression  framework  using  White's  (1982)  information  matrix 
(IM)  test  principle. 

The  plan  of  the  paper  is  as  follows.  In  section  2,  we  specify  the  model  and  derive  an 
algebraic  structure  of  the  IM  test.  Comparing  with  Hall  (19S7)  and  Bera  and  Lee  (1992) 
who  applied  the  IM  test  to  linear  regression  model  without  and  with  serial  correlation 
respectively,  in  the  ARCH  framework,  the  covariance  matrix  of  the  indicator  vector  no 
longer  has  a  block  diagonal  structure  due  to  the  inclusion  of  the  ARCH  coefficients  in  the 
parameter  vector.  The  algebraic  structure  of  the  test  is  much  more  complicated.  First  we 
derive  a  joint  test  and  then  concentrate  on  the  components  corresponding  to  the  ARCH 
coefficients.  The  test  turns  out  to  be  a  test  for  time  varying  fourth  moment,  i.e  a  test 
for  heterokurtosis.  The  test  statistic  can  be  computed  by  running  a  simple  regression, 
and  it  can  be  given  Chesher  (19S4)'s  interpretation  of  Lagrange  multiplier  (LM)  test  for 
parameter  heterogeneity.  Local  power  of  the  test  is  also  considered  here.  In  section  3,  the 
results  are  interpreted.  An  alternative  form  of  the  IM  test  is  also  computed  by  using  the 
double- length  regression  proposed  by  Davidson  and  MacKinnon  (1992)  and  is  presented  in 


section  4.  Section  5  discusses  an  empirical  example  to  illustrate  the  usefulness  of  our  test. 
A  concluding  summary  is  given  in  section  5. 

2.  MODEL  AND  TEST 


We  consider  the  linear  regression  model 

yt  =  x't(3  +  ef  (1) 

where  yt  is  the  t-th  observation  on  endogenous  variable,  x(  is  a.  k  x  1  vector  of  exogenous 
variables,  and  e$'s  are  assumed  to  follow  an  ARCH  process.  As  specified  in  Engle  (1982), 
an  ARCH(p)  process  conditional  on  the  information  set  <&t-\  is  described  as 

€t  |  $,_!  ~  N{0,ht) 

where 

ht  =  V(et  |  $,-i )  =  o-u  +  <x\£2t-\  +  a'-^?-2  +  •  •  •  +  ap£2t-,>,  (2) 


an 


d  a0  >  0.  at  >  0,  J^-i  a*  <  1-   Let  lL  =  (€?-n€?-2 e?-/))''  and  a  =  (ai5-,--  ,ap)'- 

Assuming  that  et  is  given,  the  loglikelihood  function  is  the  sum  of  the  conditional  normal 
loglikelihoods  function  corresponding  to  (1)  and  (2).  For  our  ARCH  case,  all  assumptions 
mentioned  in  White  (1987)  are  satisfied  and  the  IM  test  can  be  applied  to  this  model. 

Let  L{6)  be  the  average  loglikelihood  function  and  //  be  the  log  density  function  of 
the  t-th  observation  and  T  be  the  sample  size.  Then 

i  i 

it[0)  =  __Iog7r_._iog/lt_  _£<_ 


where  9  =  (/?',  cio.  ci'  )'  is  a  ( k  +  p  +  1 )  x  1  vector  of  parameters. 

Suppose  9  be  the  maximum  likelihood  estimator  (MLE)  of  9.    Then,  the  IM  test  is 
based  on  the  distinct  elements  of  the  matrix 

C(e)-  i^-hm     me)  oue) ,, 
c  «"  ~  f  ^Hedf  +  ( ~W )( ~W >  l 

t  =  \ 

=  A(9)  +  B(0). 


where 


«»-&%&) 


t  ^y  dOOO' 


dlt(9)    dlt(9),n 


T  ^     09  09      J 

Since  C{9)  is  symmetric.  IM  test  just  depends  on 

d(6)  =  vechC{9)  =  vech(A{9)  +  B{9)), 

which  is  a  m  x  1  vector  where  /??  =  — j^ — .  Subject  to  certain  regularity  conditions, 

it  can  be  shown  that  the  asymptotic  covariance  of  d(9)  can  be  consistently  estimated  by 
[see  White  (19S2.  p.  11)] 

i     T 

Cov(d)  =  v(9)  =  -Y,<>t(d)«dd)' 

t=\ 
where 

a,(9)  =  d([9)-  \7cl(0)A(§r1  y/,(0) 
The  IM  test  can  be  written  as 


TIM  =  Td'(0)V(0)-ld(9) 


Therefore,  under  null  hypothesis  that  the  model  (l)-(2)  is  a  correct  specification,  the  IM 
test  statistic  is  asymptotically  distributed  as  \2  with  m  degrees  of  freedom.  We  define  d(B) 
explicitly  as  follows: 

d(d)  =  [di  ,d2,d3  ,c/4  ,  c/5  ,d6  ) 


where  d\  is  a  -^j — -  x  1  vector;  dz  is  a  scalar;  e/3  is  a  ^^ — -  x  1  vector;  c/4  is  a  k  x  1  vector; 
a?5  is  a  pk  x  1  vector;  d<$  is  a  p  x  1  vector.  The  typical  elements  of  <it,j=i,...,6  are  given  below 
[  for  derivation,  see  Appendices  A  and  B]: 

1     T        1  3      l>  P  .,  P 

t-1  h=i  h  =  l  h=l 

P  P 

-{*t(€t    ~  !)  -2et)2^ah€t-ixtixt-hj  -i€((et    -1)  -2et)2-fQtfcc«-*a;*-fcia:tiJi<j,i,i«i,2,...,* 

T 


^4Bi(k-3ef  +  |)] 


,=,   2/., 


;  ^    V        T/       "*4  o-*2     1     "\~*2     ~*2     1 


t-\ 


T 


*  ■  -  f  i^W  - r<2  +  hfl *htUxUi 


+  4-(67(ef -l)-2en<]i=i,2 * 

2/i, 


T  p 

?  1  W    /     -*-i       o-*2   ,    ^  \  \~^  ~    2*        *       2*2 

t—l""  "     /i  =  l 

/~*/~*2         i\        o~*\     *  2 * 2 

-  ^(et(€<    -1)-  2ef  )*«€,_; 

+  ^(^*2   ~  1)^-j-^-jJ'=1.'2 fc,i=l,2 p 


/=1 


!/>, 


-*  ^t  ~*  €l  —  i  *  %ij  *  ^'t  —  ij 

et    —        rz-">        et-i   ~       75=?        xtj   —       7F=i        Xt-ij 


where 


ht  Vht  \/ht  y/hf 

The  variance  matrix  of  d{0)  is  not  block  diagonal.  [For  detailed  derivation  of  Var(d{9)), 
see  Appendix  B].  Therefore,  the  derived  IM  test  statistics  can  not  be  written  as  the  sum 
of  quadratic  forms  as  in  Hall  (1987)  and  Bera  and  Lee  (1992). 

As  is  well  known,  the  IM  formulaton  tests  for  the  full  specification  of  the  model  (l)-(2). 
If  we  are  interested  only  in  the  specification  of  ARCH  part,  we  may  concentrate  in  those 
particular  components  corresponding  to  the  ARCH  parameters.  These  components  are 
based  on  the  indicator  vectors  (h(6),dz{0)  and  d^{9).  One  of  these  three  indicator  vectors, 
d^{9)  is  related  solely  to  the  ARCH  parameter  vector  a,  and  now  we  formulate  a  test  based 
on  this  component.  The  other  two  components  will  be  discussed  in  the  next  section. 

To  find  the  asymptotic  variance  of  r/:J,  we  use  a  result  from  Pierce  (1981).  Let 
j/i,  t/2,  •••5?/n  be  a  sequence  of  random  variables  whose  joint  distribution  depends  on  a  pa- 
rameter 9  .  Let  9„  =  9„  { iji ,  y2, ...,  yn)  be  an  asymptotically  normal  and  efficient  sequence  of 
estimators.  It  is  desired  to  find  the  limiting  distribution  of  a  statistic  Tn  =  Tn{y\ ,  ij2, ...,  yn\ 

9n),  where  at  the  true  9  the  corresponding  sequence  T„  =  Tn(y\.ij2 Un',6)  has  a  known 

limiting  normal  distribution.  We  assume  that  for  every  9  there  is  joint  convergence  in  law 
to  normalitv: 

V  1  1        Vi  2 


^T"     )  „  v 


0 

V    V'  2 1        V'22 


and  that  there  is  a  matrix  B,  possibly  depending  continuously  on  9.  such  that 

s/Kfn  =  sfTiT,,  +  Byfii{9-d)  +  op{l) 

where  B  =  Mm  E  ($§£■).  Under  the  above  assumptions, 

iifn  ~  -V(0,I  n  -  BV2,B'). 
6 


Using  this  result,  an  estimate  of  the  variance  of  c/3  has  the  following  form: 


y  6     y~\  -*2     -*2     \/ 2*2     ~*2     \ 

V3   =   —}_^(^-^(-j)(^-let-j)t<J^,J=l 


t=l 


[The  derivation  is  in  Appendix  C:  Part  I].  Therefore,  the  test  statistic  can  be  written  as 


r3  =  <r3V3  <k 


which  asymptotically  follows  a  chi-square  distribution  with  ^^ — -  degrees  of  freedom. 

The  test  T3  can  be  perfomed  as  TR-  of  the  ordinary  least  squares  regression  of  v  on 
5,  i.e., 

r3  -  Ir'SiS'sr's'v 

where  v  =  (t)i,U2, ...,  tV)'  is  a  T  x  1  vector  with  Vt  =  {^e*4  —  3e*2  +  |)  and 


P*2     p*2 


/«T! 


5  = 


ici-i 


\^2 


e%*2 


r-icr-i 


~*2      ~*2 
6[-l  61-/) 


eT-leT-p 


1*2 


*\i  7-*l 


-p+l^l-p    \ 


^T-p+l^T-p' 


is  a  T  x  p  /J9       matrix.  [The  detailed  derivation  is  given  in  Appendix  C:  Part  II], 


3.  INTERPRETATION  OF  THE  INFORMATION  MATRIX  TEST 


From  the  last  section,  each  of  the  six  components  of  the  indicator  vector  contains  the 
common  element  vt  =  ^(^e*4  —  3£*2  -f  |).  Under  normality  and  correct  specification  of 
ARCH  model,  Evt  =  0.  .  So  the  sample  moment  rpYlt^i  Vt  would  De  expected  to  be  close 
to  zero.  Hence  a  test  for  model  being  correct  can  be  based  on  ^Yl(=\v^  a  measure  of 
sample  kurtosis. 


Here  we  are  interested  in  the  special  components  related  to  the  ARCH  parameter  a. 
As  we  mentioned  in  the  last  section,  there  are  three  such  components,  dz,d^,dQ  which  are 
related  to  the  ARCH  parameter  a  and  have  a  special  form  Xw=i(u*#<)  where  gt  is  some 
function.  Therefore,  these  can  be  considered  different  tests  for  heterokurtosis  and  each 
emphasize  the  effects  of  different  aspects  of  heterokurtosis. 

The  component  r/3  is  related  to  the  parameter  vector  of  a.  It  is  clear  from  the  ex- 
pression C{6)  that  c/3  is  based  on  the  two  estimates  of  variance  of  a.  Taking  E(vt)  as  a 
measure  of  kurtosis.  d;i  measures  the  change  in  the  kurtosis.  more  precisely  it  tests  whether 
the  kurtosis  depends  on  the  cross  product  of  the  lagged  residual  squares.  Following  Chesher 
(1984),  we  can  also  give  a  test  for  heterogeneity  interpretation  to  c/3.  Suppose  the  ARCH 
parameters  a  are  varying  around  a  mean  with  finite  variance.  This  can  be  formulated  as 
at  ~  (cv,Q).  Then  X3  is  the  LM  test  for  testing  Hq  :  £1  =  0,  i.e.  it  tests  the  randomness  of 
the  ARCH  parameters. 

Next,  c/5  is  based  on  the  relationship  between  l3  and  d.  The  assumption  of  a  symmetric 
distribution  implies  Ee*1  =  0,  Ee*  =  0  and  that  allows  us  to  omit  the  second  part  of  j5.  The 
third  part  can  also  be  omitted  because  it  deals  with  heteroskedasticity  arising  from  the  cross 


products  between  e*      and  x*_  -  and  is  not  very  important  from  a  practical  point  of  view. 

Then,  d5  reduces  to  d5  =  ±  Y,T=i  Kl^*4-3^*2"!  )  Sf=i  &l^*-lxt-liit-j)i=i,2,...,k,j=i,2,...,p 
and  this  describes  a  relationship  between  ut  and  ]T)jLj  a(e*_/^*_^^*ij  and  allows  us  to  test 
heterokurtosis  caused  by  cross  products  between  lagged  error  terms  and  lagged  exogenous 
variables.  d§  is  an  expression  arising  from  the  two  estimators  of  the  covariance  between 
<72  and  d,  and  can  be  used  to  test  the  heterokurtosis  due  to  the  square  lagged  errors. 

Since  the  covariance  matrix  V  is  not  block  diagonal,  these  three  test  statistics  are 
definitely  correlated  with  each  other  and  also  with  the  other  three  components.  To  get 
overall  test  of  the  model,  it  is  necessary  to  have  a  joint  specification  test.  This  can  be 
obtained  by  using  the  results  in  Appendix  B. 

4.DOUBLE  LENGTH  REGRESSION  FORM  OF  THE  TEST 


Davidson  and  MacKinnon  (1992)  proposed  a  double-length  regression  (DLR)  to  per- 
form the  IM  tests  on  models  which  can  be  expressed  as 

ft(yt,0)  =  et,      ,=  i r,      et  ~JV(0,D- 

For  this  class  of  models,  the  contribution  to  the  loglikelihood  function  from  observation  t 


is 


It  =  --Io(j{2tt)  -  -ff  +  kt, 

where  kt  =  log  |  |^-  |=  log  |  /]  |  is  the  Jacobian  contribution  to  lt.  The  DLR  uses  2T 
"observations."  The  regressand  is  ft  for  "observation"  t  and  one  for  "observation"  T  +  n, 
n  =  1,2,...,X,  and  the  corresponding  regressors  are  respectively  —  Fti  and  Ktii  where 
Fti  =  -A:  and  Kti  =  jnfr-  The  test  statistic  is  then  the  explained  sum  of  squares  from  this 
artificial  regression.  In  order  to  obtain  the  DLR  form  of  the  IM  test,  an  explicit  alternative 


hypothesis  of  the  model  ft{yt,Q)  =  £t  i-s  needed.  Chesher's  (1984)  result  suggests  the  use 
of  the  following  model: 

/*(V*,0  +  Ct)  =  c<,€i  ~  iV(0, 1),  (3) 

where  #  =  (,^',Q'o,a')'  is  a  k  +  1  +  p  dimensional  vector  ,  and  Q  =  (0',0,a/{)'  and  ??* 
are  k  +  1  +  p  and  p  dimensional  random  vectors  respectively,  with  i]t  being  distributed 
independently  of  et  and  of  rj3,  s  /  t.  //<  ~  iV(0,  2S7).  By  taking  a  second  order  Taylor  series 
expansion  of  (3)  in  Q,  we  have 

ft{y,.0)  =  e,  -  F,Ci  -  7}CF*0  =  e'  ~  Filt  -  ^ItFtlt  =  (Jt  {say), 

where  Ft  is  a  1  x  ( A"  +  1  +  p)  row  vector  with  typical  element  -^  and  Ft  is  a  (k  +  1  + 
p)  X  (k  +  1  +p)  matrix  with  typical  element  de  [^  ■  Ft  is  a  1  x  p  row  vector  with  typical 
element  ^-  and  F*  is  a  p  x  p  matrix  with  typical  element  da  ^Q.  ■  Note  that 

and 

VW^  |y<)  =  l  +  2tr(ft-Fl7,F,), 

where  UT"  denotes  transpose  of  a  matrix.  Thus  locally  in  the  neighborhood  of  Q  =  0.  the 
model  is  equivalent  to 

ftiyuO)  -  E((j,  1  tjt) 

=  i)t,   0t  ~  N{0, 1).  (4) 

i.e, 

/)(y,i)  +  />ii)f;)      , 

qtiyt.tiAl)  =  r  =  Ut. 

(l+2tr(QFtTFt))i 

The  loglikelihood  function  of  this  model  can  be  written  as 

1  ,       ,  1    2 

/,  =  --Io(j(2tt)  -  -qj  +  rt, 

10 


where  rt  =  log  |  |^l  |.    Now  a  LM  test  for  Q 
from  the  following  regression 

//i\       /Quj\ 


0  in  the  above  model  can  be  computed 


1 


Qtn 


{QtJ)  +LO< 


(Clij)+wt 


w 


V  1  /        \RTijJ 

here  i  =  l,...,p,j  <  /',  and  £t  =  (e'/_j , ...,  €j_p+1e^_p)  is  a  1  x  /;(/;^"      vector.    [Detailed 
derivation  is  given  in  Appendix  D:  Part  I  ]. 

Since  the  loglikelihood  functions  of  the  Davidson  and  MacKinnon's  and  Chesher's 
models  are  the  same  under  Hq  :  Q.  =  0  although  they  differ  under  H\  :  Q.  =£  0,  the 
constrained  MLE  obtained  from  those  loglikelhood  functions  are  the  same.  And  the  score 
vectors  of  those  loglikelihood  functions  at  the  constrained  MLE  are  also  equal.  Using  the 
terminology  of  Godfrey  (1988,  p. 72)  we  can  then  say  that  the  models  (3)  and  (4)  are 
locally  equivalent  alternative  (LEA)  to  the  ARCH  model  given  in  (l)-(2).  Hence  the  LM 
test  statistic  for  Ho  should  be  the  same  for  both  the  loglikelihood  functions. 

Given  the  above  results,  the  theoretical  local  power  of  T3  test  is  equivalent  to  that 
of  the  LM  test  with  the  Davidson  and  MacKinnon's  model.  To  derive  the  noncentrality 
parameter  of  the  LM  test  we  consider  the  simple  case  of  p  =  1.  Note  that  now  we  have 
0  =  {j3',a0lai)'  is  a  {k  +  1  +  1)  x  1  vector,  ;"  =  c*i  +  r)  and  rj  ~  JY(0,2<72).  For  the  LM 
test,  the  null  and  alternative  hypotheses  are 

H0  :V(,))=2a2  =0 


11 


and 

6 


Ha  :V(n)  =  2cr2  = 


v/T 

and  similarly  for  the  IM  test 

Ho  :  h,  =  a0  +  (-\ie2_l 

and 

Ha  :  ht  =  Qu  +  (ai  +  ij)e2t_l,  (5) 

where  rj  ~  iV(0,2<72).  Thus,  locally  in  the  neighborhood  of  a2  —  0,  the  model  is  eqivalent 


to 


,        fl(y.0)+y2Ft:y, 

<lt{yt.0*<r~)  =  ,  „    „ — r  =  /'■'■ 

(1+2^F«F«)» 


with  log-density  function 


/,  =  --log(2w)  -  -q2  +  rt. 


The  variance  of  the  test  can  be  obtained  using  E[—  .-^J2'2].  Therefore,  the  test  statistic 
under  the  local  alternative  follows  a  noncentral  \2  with  1  degree  of  freedom  and  noncentral 
parameter  82 ^E^2t=[  —  -^-Ai  ■  [Detail  derivation  is  in  the  Appendix  D:  Part  II.]  The 
variance  we  got  here  is  the  same  as  we  obtained  for  T$  in  Appendix  C. 

Following  the  procedure  of  Engle  (19S2).  it  is  easy  to  find  the  second  and  fourth  mo- 
ments of  a  first-order  random  coefficient  ARCH  process  i.e.  of  model  (5).  These  moments 
will  give  some  idea  about  the  nature  of  the  implicit  alternatives  for  the  IM  test.  Letting 
A,  =  (e?,e2)', 


The  condition  for  finite  unconditional  variance  is  same  as  in  the  non-random  coefficient 
case,  that  is,  0:1  <  1.  while  to  have  a  finite  fourth  moment  it  is  now  required  that  3(a2  -f- 
2<t2)  <  1.  For  the  standard  ARCH(l)  model,  this  condition  is  3oj  <  1.  If  these  conditions 
are  met,  the  moments  can  be  computed  as  [For  derivation,  see  Appendix  D:  Part  III]. 


E( 


4\_r       3a0       n  1-a'f 

et  1  ~ 


1  -Ql)2Ul-3(a2+2<72) 
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£(€7)    = 


1-or, 

Therefore,  the  kurtosis  of  the  random  coefficient  ARCH(l)  model  is  3(1  —  a\)/[\  —  3(a2  + 
2a2)]  which  is  higher  than  3(1  —  of ) / [  1  —  3oj],  the  kurtosis  for  the  standard  ARCH(l) 
model,  for  a2  /  0.  Since  the  unconditional  variance  remains  the  same,  our  alternative 
hypothesis  can  take  acount  of  higher  degree  of  nonnormality. 

5.  AN  EMPIRICAL  ILLUSTRATION 


Engle(19S3)  estimated  the  following  reduced  form  equation  for  inflation  using  quar- 
terly data  from  1947-IV-1979-IV 

Pt  =  fa  +  MPt-i )  +  MPt-2  +  MPMt-i ) 

e,  |  $(_i,~  -V(0,/?., ) 

i=i 
where  P,PM,W,  M  and  t  are  respectively  the  rates  of  change  of  the  GNP  deflator,  the 
rate  of  change  of  the  import  deflator,  the  rate  of  change  of  wages,  the  rate  of  change  of 
money  supply  and  a  time  trend.  This  inflation  equation  includes  two  lagged  dependent 
variables  and  the  conditional  variance  h(  is  assumed  to  exhibit  a  two-parameter  ARCH 
process  of  eight-order  with  linearly  declining  weights. 

As  reported  in  Engle  (19S3),  the  above  model  satisfies  the  standard  diagnostic  checks. 
Also  Bera,  Higgins  and  Lee  (1990)  found  that  above  model  passes  the  Pagan  and  Sabau 
(1987)  test  for  correctly  specified  conditional  variance.  The  T3  test  examines  this  model 
at  p  =  1.2.3.4.5.    For  p  =  1  and  2.  the  values  of  T:i  are  respectively  0.445  and  6.569. 
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These  are  not  significant  at  both  1%  and  5%  significant  levels  for  which  the  asymptotic 
critical  values  are  6.635  and  3.S41  with  1  degree  of  freedom,  and  11.345  and  7. 815  with  3 
degrees  of  freedom  respectively.  However,  for  p  =  3,4,5,  the  test  statistics,  61.98,  139.60, 
and  199.34  respectively,  are  highly  significant  at  both  1%  and  5%  levels.  This  may  imply 
that  T3  can  reveal  model  misspecification.  while  the  standard  diagnostic  checks  may  fail 
to  do  so. 

5.  CONCLUSION 


In  this  paper,  we  have  provided  an  application  of  IM  test  to  the  linear  model  with 
ARCH  process.  We  give  the  computation  and  interpretation  of  the  resulting  test.  Because 
ARCH  error  are  involved,  the  framework  of  the  information  matrix  test  is  much  more 
complicated  than  those  derived  by  Hall  ( 1987)  and  Bera  and  Lee  (1992)  for  simpler  models. 
In  our  case  the  variance  of  the  indicator  vector  is  no  longer  a  block  diagonal  matrix,  and 
therefore,  the  components  of  the  indicate  vector  are  not  asymptotically  independent.  We 
applied  one  component  of  IM  test  to  the  Engle  (19S3)  model.  The  test  we  use  above  may 
have  higher  power  than  the  tests  which  Engle  (19S3)  had  used  for  his  model  in  the  sense 
that  the  proposed  test  rejects  the  model  specification  while  other  diagnostic  checks  can 
not. 
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APPENDIX  A 


The  Derivatives  of  the  Log-likelihood  Function 

For  our  model,  the  parameter  vector  is  9  =  {[3' ,  cto,  a')'  and  the  log-likelihood  function 
for  the  t-th  observation  conditional  on  the  information  set  $t-i  is  given  by 

1(0)  =  ^0</27r  -  -loght  -  T^j-t]- 


The  first  derivatives  are 


da        2h,    h, 
and  second  order  derivatives  are 


dU         e]     1    dh%         1    dht       et 
dj§  ~  ^h~t2ht  d(3  ~  2ht~d3  +  T*' 

dh  1     c? 

— -  =  — (—  -  1) 
ctao       2/?t    /?/ 

91,         1  ,«?      ., 


<92/,  e<      <9/?,        e<  <9//,    ,        1        ,        1    t2t   dht  dht 

Olid?  ''  ~~k\XtW'  ~  Tti~oiF'1'  "  h~tXtXt  "  T^WtWW 

1      e?  0/?,d/?,        e?       ,%   1     <92/i, 


2/Y2%         '<9/i  dp        It,        ''2htdSdf3 


d2h  I   A      „        1    c? 


da0<9rvo  2/7?    /*,  2/??/?./ 

1      ej  ,1 

-—(—-1)6,^     ~  r—;—€tei 

2hj    ht  2hj  h, 


d2h  1    ,  e?       _       ,        1    67        , 


dad 


d2lt  1      e?  dht         I    e^dlh        1_ 

djdao  ~  ~2hi{Yt~l)!hJ    '2iiih~t!)J    'l^et'Vt 

d2U  1    ,e?       ..a^t     f         1    e?3/*<     , 


d£da'  2/if v  ht        'dp-      2/Y2  ht  dj3~ 

1  1  ,e2  d2/i, 


/if        -       2/i,/i,         'oWo' 


0ao0a'  2hyht  -       2hjh, 

1G 


w 


here 


and 


The  outer  products  are 


-  =  -  ^latet-iXt-i 


03 


i=i 


d'2ht        A 

7  =  2^  2a«x<-«x*-« 


d/Jd/i 


i=i 


a/t  a/, ' 


2       1       C)/?/   5/2  <' 


a/?  d/J 


fcl 


4/??  #d  a/? 


+ 


2/i?  v  /i 


f--D^: 


<  '«/ 


a/r* 


2/i 

cM,  dlt  '         1 


dp1  da0        4/i 


1 


dhdh' 

d/3  do~        4fc 

d/«    dZ,  '         1 


<9a0  5oo         4/? 


1 


dlt_dl^  _ 

dcio  5a  4// 

0U_dW  \_ 

da  da  4/t 


q  _      20/i, 


1 


h 


0 


0  2/?7     A  # 


V 


—  -  l)2  — 


t-"2 


P--DV 


'  i  \2  ' 

7 1)  e«c<  • 


+^t4t'i)etx^ 


17 


APPENDIX  B 


Covariance  Matrix  for  the  Information  Matrix  Test 


A  consistent  covariance  estimator  for  the  IM  test  proposed  by  White  (1982)  is  stated 
as 

V{0)=  ^22at{9)at(0)'  (B.l) 

where  at{9)  =  d(9)  -  \/d\9)A{9)~1  y  U{9).  Let  us  begin  with  the  indicator  vector  d{9) 
which  is  defined  as 

d{9)  =  vech[C{9)]  =  vech[A(9)  +  B(9)] 

where 

i,«,-lf[?Mi     -if 

A^>-  ]-LI  d0Q0,  !»=«-  rZ. 

l=\  t=l 

2et        cl    I     2«f  c    -,./  ■*(■*>  2  i      1    Jj  /*  1    ,.    c  (t^t  1    .-    C  c  I         ttxt  r  I  1   f.    ,, 

-£f  x^<  +  *«**«  -  —  ~  /TJK'\/  +  7^/0     jtM^i  -  iTT-      hj^tStei  -  -fy-et   -  j-<V7* 


and 


(  t  '  !  I 


7'  T 


1=8 


t  <  *  t 

where 

9  > 

/?.,        // , 

IS 


1=6 


9 


ht 

v 
y    otjEt-jXt-i  =  Zt 

i=\ 

P 

Y^otiXt-ix't-i  =Q 

1=1 

p  p     p 

i=i  ,=  i  >=i 

and  //;  is  a  k  x  p  matrix  with  (i.j)th  element  €t-t.vt-,j ;,i=i,2 />  and  7=1,2 k-  From  A(B) 

and  5(0),  C(0)  is  derived  as 

1     I 

C(0)  =  -4(0)  +  5(0)  =  -^ 

<-i 

1   \  P'        l_ii  >•'  i    \  i    \  c  i 


where 


2€«       e/  ,  e? 


[.r,.r',  +  ^a,.r,_,.r/,_,]  =  /t, 


Therefore.  c/(0)  is  given  by 


r 

/=i 

where 

(I, (6)  =  ( </', , .  d'n .  rfj3 .  d'iA .  r/',5 .  (/'„; )' 
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d(B)  =  jJ2dt(d)  ( 


and 


dt\ 


=  {(£?  -  1)[*J?  +  J2  *i*£ii]  +  2(  Je?  "  3e?  +  2 


j=i 


2^  «i  Ct-i^i-il   +  2  2^  2^  <**<*?  ^-i^-j^-il  Xt-jl  ) 

P  P 


i=  i 


»=i 


(^- 1)1*3 +  £>*&*] +  2(Js4-3ef»  +  5 


1=] 


(X  *   '^   ^* '^  *'^  »"\      x  x  rt         ~         ^  *  ^  :*:  *  *  \ 

2^Qiet-txt-ik  +  2LLa'aJe'-'e'-rr'-^-i''-^-) 

i=l  Kj   7=1 

p  p 

+  2(2eJ  -  e*(et*2  -  1))(^  ajej.^ar^.  +  ^  aje*_-a;t*_tifcxt*fc), 


»'=] 


(=1 


d'2  -  i)(.«-;,.«72  +  £>*?_,-,*:_»]  +  2(^;'  -  3ef  +  | 


i=l 


(2^ ai  €t-ixt-axt-i2  +  2  2L  2^  n'a^'-'e'-vr'-'i-r'-j-2 ) 

/'  /' 

+  2(2eJ  -  e*^2  -  1))(]T  afi-ix*ixt-i2  +  ^Z^*-^*-*!1")'-' 

*=i  i=i 

(^'2  ~  lK(fc-i)^t  +  E^'-r'-'(A-i)J*-iA-]  +  2(ir4  -  3^2  +  |) 

/=i 
p  /' 

{2_sQiet-iXt-i{k-l)Xt-ik  +  2  2^  l^Qi0iJet-i€t-jXt-i(k-l)Xt-jk) 

i=l  i<j    ;=1 

P  l> 

+  2(26?  -€*t{t*2  -  l))(^0|Ct-,-.rJ(ib-i)a?r-,-fc  +  ]T  Qr,-€7_iarJ__i(ib_1)a:JJfc)]' 


(=i 


/=! 


is  a  -^-tt — -  x  1  vector. 


dti 


1    .1 


*l  o;*2 


-e,     - 


2/,.?v2  ' 
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3?'  +  o) 


is  a  scalar, 

5(5^-3^  +  1^^ 

is  a  i^±^  X  1  vector. 


is  a  fc  x  1  vector. 


11  q      p  1 

<*t5  =  ["(^e<     -3et    +  -)2_^al€t_rvt_letLl-  ~(2et  -et(et-  -l))xtetll 


1=] 


1 ,  -.,2 


--(e,     -l).r,  _,€,_! 

11  3    J-  1 

/       ~*1         o  2*2     1  \  \    ^         ~*  *        ~*2  /  •■» ~*  ~*  /  ~*2         1  \\     *  ~*2 

--(-ef     -36,     +  -)2^ai€t_rrt_pet_p  -  -(2e,  -  et(et     -l))xtet_j 


t=i 


/  ~*2  1  \     *        ~*       V 

--(et    -  l)aj|_,Cf_p] 

is  a  kp  x  1  vector, 

11  ^  11  3 

o«6  =  lrf(o€«    ~3e     +o)6'-i T7~(oe<    ~3e'    +o)6'-/>! 

2/?/    -  -  2ht    -  - 

is  a  p  x  1  vector.  Next  we  consider 

Vd(e„)  =  nhm-^£;[-^-]. 

Using  the  normality  assumption  of  the  e(  and  taking  expectation  conditional  on  the  in- 
formation set  $/_i  iteratively,  after  some  algebra  we  can  get  the  following  simple  form  of 
W(#o) 
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Vd(0o) 


where  Xjdn  —  (m.x'n, ...,m.TjtJt,?Ti.Ti2. 


/      0  V^i2  V^13\ 

0            0  0 

0            0  0 

0            0  0 

V4n         0  0 

V    o        o  o    / 


n?.r(^._! )/,.)'  is  a  -^ — -  x  1  vector  with 


T 
1  A    1 


mxi-i  = 


i<j.i,j=l,2 k 


i  h'i 


t=i 


i    :.    „     *■(*■+!) 


V^i3  =  (wii,...,iyjtjk,it'i2 u'u— i)*)'  is  a  — -, —  x  p  matrix  with 

1     -     1  £ 


i=i 


i</,t,j=l,2 A- 


and  V^5i  =  (~ii-  -22 -/>/>)'  i-s  a  M'  x  ^'  matrix  with 


7' 


•r^c  T  ^  77[f'-«J,'-'(Z  aJc'-J*l-j )] 


T  tl  « 


7=1 


1=1.2 /> 


This  implies  that  Sjd{0u)  can  be  estimated  consistently  by  the  \jd{9)  which  is 


Vd(*)  = 


(  ° 

W12 

Wn\ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Wsi 

0 

0 

V   o 

0 

o    ) 

where  for  example,  V^r~>   =  (rnxu mxkk*mzi2 '"^'u-i  )k  )'  is  a 

with 


k(k+i) 


x  1  vector 


/??  x  i 


1       T  '' 

1  "<    f=l  /=l 


09 


0 

0 

1 

2h} 

--*-e' 

2/it2-< 

1     € 
2/»2-< 

--±-e  e' 

i<i,i,j'=l,2 A;- 

Similarly,  we  can  simplify  .4(0)  as  follows: 

A{6)  =  |  ° 

0 
Finally,  from  Appendix  A  S/lt{0)  =  §#  is  given  by 

/-(«*!-i)ELi«i^;-,+^; 

Given  the  above  expressions  for  d{8),  S7(l{9)<  -4(0),  and  S/l({0),  using  the  formula  (B.l) 
we  can  obtain  an  estimate  of  the  covariance  matrix  for  the  IM  test.  Unlike  the  cases  of 
Hall  (1987)  and  Bera  and  Lee  (1992),  here  .4(0)  is  not  block-diagonal,  and  this  results  in 
V{6)  to  be  non-block-diagonal.  And  the  final  expression  of  V{6)  is  very  complicated  and 
is  omitted. 
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APPENDIX  C 


Part  I. 

Here 


Using 


we  obtain 


Defining 


We  have 


i<j,i,j=l,...p 


V(a:)  =  £[r(j:|y)]  +  V[f7(a;|y)]1 


r 

V(</.,)  =  vi,  =  ^[^^;:(^)'(6;:(^)G]. 


1,1   =   r  l^(et'-'€,'-J)  ^t-i^t-j)^ 


Also  note  that 

+  (2  '  '    +2,/"1     '~J    df3         ht[      '  '^   *-*-J 


Ir1^         o*2         3    d/li 


*2      *2 


^3  *   /         * 4     i    o    *2  \    *2      *2  -  /  ^     *  I         o    *2     i     "  \    *2    „*2 

7J—  =  7~(-^    +3e,    )e,_,-e,     ■-— (-e,     -3e,    + -)et_let_J 

OQq  llt  lit     -  2 

^"3  1    /         *4     i    o    *2  \         *2       *2  -   /  ^     *4         o    *2     ,     "  \        *2       *2 

doT"/^"6'    +3e,-)£i€<l1e/ij-—  (-e<     -3e<    +  -  ^e^e^ 
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Op  oa  o  da 

Part  II. 

The  test  T3  which  we  got  in  section  2  can  be  written  as  nR2  form  by  running  a 

regression  0  on  S. 

1 


ch  =  ^S't), 


where 


/.'\ 

O9 


\  VT  I 

is  a  T  x  1  vector,  vt  =  (  ^e*4  —  3e*2  +  |)  and 


5  = 


~*2      ~*2 
6l-lel-/> 

e2-le2-p 


2*2  2*2       v 

tl-p+ltl-p    \ 
2*2  2*2 

e2-p+le2-/> 


\   2*2      2*2  2*2      2*2  2*2  2*2        / 

veT-l€t-l       •"       eT-\et-i>       •••       eT-},+  l€T-p/ 

is  a  T  x  p(prfl)  matrix.  From  the  result  above,  we  know  that 


V(d,)  =  E(V(d,  \*t-iY. 
1 


AT 
0 


4T2 


rS'E(tw')S 
S'S. 


T3=d'3V(d3r1d3 

0 


TR2  =  T 
1 


PSiS'S^S'v 


v'SiS'sr's'v 
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v'v 


plimT—oo-=-  =  6 
T3  =  TR2, 

where  R2  is  uncentered  coefficient  of  determination  of  regression  v  on  S. 
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APPENDIX  D 

Part  I. 

To  derive  the  double  length  regression  (DLR)  form  of  the  IM  test,  we  first  note  that 

IJt   —  X'fii  =  €t,         ft    |   $/-!    ~  iV(0,  /?f  ). 


and 


ft(yt,0)=  V'    J^j  =  /*,.     ^  |$/-i  ~iV(0,l). 
s/hi 


Following  Davidson  and  MacKinnon  (1992),  a  locally  equivalent  model,  under  the  alter- 
native hypothesis  that  Var(a)  =  2Q.  can  be  written  as 

"  fl0.        f,(y,,e)+tr(nF(*) 

q,{yt<V,U)  =  r— r  =  0f. 

(I  +  2tr(nF?  Ft))? 

Using  the  notation  as  defined  in  section  3,  we  can  write 


Ftl 

Of, 

h    .2 

„  a  €t-r 

2/i? 

Ft,  j  = 

d2ft 

da,daj   le=e 

3e,  „2     „2 

„    5    €t-l€t-J 

?,      0ft 
h  -  On,  ~ 

1 

vC 

Ku  ■ 

dlogf't 

Da,      lo=«~ 

1      .2 

2/i? 

Ktij  = 

d2locjf't 
dcxida-j  lo=& 
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3    .,    .2 

^   5   €t-l€l-j 

4/i? 

Qtij   —  Qtij    \Q-Og-0 

Oqi    | 
=  Ftij  —  ftFtiFtj 

o~  ~3 

Oct    „2       ->  ct      -2       ~2 

—       „  s  et-i€t-j  ~et-i€t-j 

Hi?  4h? 

4/i?    Vht       hj 


zn*  ■'■     ,     '^€<  et    m  -4  -2       -2  -2  -2       \ 

£*  =  T7T(^^  "  rr)(^-i et-i€t-2 et-p+iet-P) 

4/'F     Vht       h] 

1    ,   3««         c?     - 


4/>2     V%        /if 
where  e<  =  (e|_!, ...,  €?_1eJ_2 ef_  +1e,_   )  is  a  ^-^ —  x  1  vector. 

dlogq', 
=  A'jjj  +  KuKfj  -  FnFtj  -  ft{K,tFtj  +  FtlK,j) 
4/*2    y/Ki       ht  s/Ji    ht  3 


13  1  '">      e2    - 

J??  =  -^-(-7=  +  ^  -(1  +  -^=)^-)^. 
4/?f    Vh,       ht  yjht    ht 

Given  the  above  results,  the  DLR  form  of  the  test  can  be  expressed  in  terms  of  the  following 
regression 


i\K*>+w' 


,/,r  ^e' 


«  Q(J+cj, 


1     /      3 


+  f  -(1+     '    )|i-)e,. 


1/(r    v/^7       /('  >A7   ''« 


2S 


Part  II. 

To  calculate  the  non-centrality  parameter  for  local  alternative,  we  consider  the  simple 
case  with  p  =  1.  The  model  under  H^\  a1  =  0,  is 

r,         m  lit  ~  Xtfi  *t 

jt{yt,v)  =  — -?==—  =  —7=  =  nt 
s/ht  v'n 

and  \it  ~  Ar(0, 1).  The  loglikelihood  of  this  model  is 

i.  ..      i..2 


and  therefore, 


lt  =  -log{2ir)  -  -ft  +  kt. 


Fn  = 

Of, 
dc\  i 

€(      2 

2hl€t-1 

Fm  = 

d2ft 

daj 

3e<    4 

5    €t-\  ' 

4ft? 

/; 

Of, 
Oy, 

1 

^/i  = 

0Fn 

Oy, 

1    2 

2/i/ 

n, .  = 

dFiU 
dyt 

-     3    e< 

—              5   tt-\ 

df, 

et 

dh,        xt 

03         2hf  03       s/17t' 

df[  1    Oh, 

03  ~      ojA  d(3' 

0Fn         3e,    .y     0ht        xt     .,  et 

93        Ah]         °P       2hf  hf 

OF;,         3     2     Oht        1 
°j         4/?.,2  ^        fc/ 

<9Fm  15f/   .,      Oh,        3.r,    4  36/   3 

D,i  Shf  di       Ah]  hf 
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dF;u            15  4     cM*        3    3 

°P             Sh-f  °P        hi 

Oft  ft    dh.t 

da  9/,f  da 

-"  t 

df[  1    dht 

da_  9l|  <9a 

<9Fn  3e<    ,     cM, 


0« 

~              5    W-l     O         • 

da 

3     .,     % 
4/i?          ^ 

dFtu 

da 

15e<    j      5/i< 
Shf          da- 

dFln 

da 

15  ^     cto, 
S/i?         ^ 

Here  a  =  (qo,Qi  )'.  The  model  under  H„ :  a-2  =  -4=,  is 

,       a     2,       //(.'/ii)+^Fni 
(l+2^FMFn)i 

and  f/j  ~  Ar(0, 1).  The  loglikelihoocl  function  of  this  model  is 


1,       ,  1    , 

/,  =  --i0(J{2n)  -  -q~  +  r, 


We  denote  the  information  matrix  of  the  model  as 


"  U-2]     v22 


where 

0(1  dp'         dp  da 

Vl1   :       _A  c)2/,  0*1, 


dadp         dada1 

is  a.  (A:  +  2)  x  (A-  +  2)  matrix, 

V12  =  K]  =  -£l      .)2ft 
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is  a  (k  +  2)  x  1  vector,  and 

d2lt 


Vio  = 


is  a  scalar.  Now 

dqt  _(l  +  2a2FnFn)i(|^  +  ^^)-(/f+a2^)(l  +  2<72FnFni)-i2(72F<1^ 
dp  ~  (l+2a2FtlFtl) 

dgt        (1  +  2<r2 FtiFti^Fm  -  (ft  +  <T2Ftn)(l  +  2<T2FtiFnr±FtlFti 
da2  (l  +  2a2FtlFn) 

dgt  .         _  dft 

dp  '*2=0~  dlV 


"7"    |<72=()   —   —  TT^'ll   +   -TTTjt^tl^tl 


dp  da2  x°  "u  dp    llx       dp 

_  3e2    .,      <9/?<       3artei  _4 


Sir}  (~[  OP        4/ij 


3  w-i 


.,    ^-1- 


S/zf   '_1  dp         Irj 

Ft    dqtdqt  1       1*1-    3  ,«    a/"       3  ,4    ^/?'-o 


/        ^m      w"< 


a2<;«  ,        9f„,  aF„  a/, 

aw  l"2=0~  ~aT  ~  -/"F" ~alT  "  F,lF"  W 


f;' aw7  l-'-°  =  ~/( ~W  +  2f> F"  ~W  +  /"F"F"  a"* 

15e2   4     dh,       3etx(   .,  3e2   3 

H — r7T-e/-i  +  -rr6/-i^t-i 


3e'J    4      d/?.*       c Jar  1  4  e|_  3 

*h*  Mdp     2/1J  '-1 "  /?;  '-1  '-1 


e*     j      a/?,       ejar(   4 

'  t   y 


shf'-]  dp      4/if  '-11 

IP/      d'2fr   I     l«    ,       3  -I   0h< 
E{~qtdpd^l<r2=ol^-i)  =  ~^e^W 
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Then 


dr<        r/i    ,   o    2rp    r    .l,dft  2dFtll\ 

—  =  [(l  +  2aFnFn),(—  +0-^-) 

OF 
+  (ft  +  *2Ftn){l  +  2<T2FtlFtlrh(T2Ftl-^j- 

OF' 

-  (ft  +  <J2 Ftn){l  +  2a2 Ft[Ft[)-±2a2  — £-Fn 

op 

dF 
-{ft+cT2Ftll)(l+2a2FnFtl)-L>2*2F'n^± 

Of,        2OF,u  2  v-i«  2  t?i  rp 

-{^r-  +cr- -——){! +2a- FnFn)    ?2cr-FnFn 
03  Op 


+  (fl+(T'Fin)(l+2a'FnFn)-i2<T*F'nFn2<T"Ftx 


03 


(l  +  2a*FtlFti) 


Q2r, 
030a1   l<72=°  ~ 


1   OFJn       0fi'       OF^      ^fV       OFn 
ft'    dfi        ~f,     n  dp       ~f,     "   03 

~ft,Jn    n  d(3       f?    tnd{3 


Shf'-1  03     '  h-f'-l'V'-l  +  4h!€'-[  0;J 

e[  3  ej  <9/>,        etart  4 

3     .,      Oh,         e2     4      Oh,       3e2    4      d/i< 

+  S^f'-'a7   'Ih~ie'-[  ~03+Ih:i€,-1~03 

etx(   .,  e2    3 


£(0^u*=ol*'-1,-4fc?e'-1aJ 


^-^fe  +  fe*'*'-0'*'-1]-0- 
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Since 


02lt  02q,  dqt  dqt  02rt 


■2    "        Cb  PiQ5\-1  ?\a    Pl~1      ' 


0,30a2  ll  Odder2       0(3  da2       0/30a2 


02lt 
OlWa2 


E(ttt^  |^«ol  *«-i)=0, 


Using  the  same  procedure,  we  can  obtain  the  E(  jrrtt  U2=ol  $t-\ )  —  0.  Then,  we  have 
V\2  =  V2'i  =  0.  Because  the  inverse  variance  of  the  non-central  parameter  is  equal  to 
V22  =  (V22  —  J^iVn  Vi2i)""l>  depending  on  the  above  results,  V22  =  (V^)-1-  Now  we 
calculate  V22,  using  the  following  derivatives 

dgt   =  (l  +  2<72FMF/1)iFm  -  (ft  +  v2Ftn)(l  +2<T2FtlFti)-iFtlFtl 
da2       '  (l  +  2a2FtlFn) 


r)2 

■777  U 2=0=  —'2Fti\FtiFti  —  Zft{Ft\Fti)~ . 


0(a2)2 


d(l<    d(l1     I  I  r  rr     r     *2 

TT~. T7T-7     <r2=0    —   l-T/1  1    ~~  Jtr  t\t  ,\  J 

do-2  C/<72 

=  F2n  -■2flFtnFnFn+ff(FnFn)2. 


,        dqt 
(lt  = 


(l+2g2FnF„)i(/;  +  <72F/11)-(/f  +a2Fm)(l+2<72FnF/1)-i2(72F;iFa 

(l4-2^F„Fn) 


?•,  =  /o</</| 

-/o^/[(l  +  2<72F/1Fn)"(/;4-a2F;il) 
-(/, +a2Fm)(l+2a2FnFnr^2a2F;iFn 
-/o^(l  +  2a2FMFM) 
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da2 


=  {(1  + 2a2  Ft[Fnrnf'l+ a2  F'tn)FtlFtl 


-I«      9 


+  {l  +  2crlFtiFn)*F'tu  -  Flu(l  +2a2FnFtl)-^2a2F'nFtl 
+  {ft  +  cr2Ftll)(l  +  2(j2FllFtl)-*FnFn2a2F'tlFtl 

-  (ft  +  a2Ftll)(l  +  2a2FnFn  r^^F,,] 
/[(l  +  2a2FtlFll)Hft+<72F'tn) 

-  (ft  +  a2Ftu)(l  +  2a2Ft]Ftl  )"ha2 F'tlFtl\ 
-2FtlFti/(l+2<r2FtlFn). 


pp.  -I 

«,    9V>  l^2=°  =  ,  £  i\o  Ui(-.t't(Ft\Fn)-  +  F'nlFtiFn 
oyer- )  [jt  )- 

+  FtlFtlF'tu  -2FinF'nFt]  +  2flFtlFtlFl1Fn 
-2FtnF'nFti  +2flFtlFllF'nFtl) 

-  (f'tFnFtl  +  F'in  -  2flF'nFtl)2}+MFnF,ly 

=  jp[-ft2(FnFn  )2  +  SflfdFn  )3F'n 
Jt 

-  +f'tFtnF'nFt[  -  F'inF;n  -  IfhF^Ftlf 
+  4ftF'luF'nFn]+4(FnFtl)2. 


02q, 


d(«2) 


7-7  Ua=o=  —~}tFt\\Ft\Ft\  +?>};{Ft\Fti 


d2l, 
0(a2)2 


d(o2)2       'da2'         0(a2)2 
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and 


d2l 
o,    2 ,2  Ia2=u  =  4ftFtuFtiFti  —4=ft(FtiFti)~  -  Ffn 


f 


'  /   Z7      \3  77' 


1       „, 


+  $^(FtiYF'tl-4-FtnFl1Ftl-—7Fl11Fi 


f 


ft 


ft 


(ii 


-  ilL(F'nFnf  +  4-^F'luF'nFn  +2(FnFnf 

ft  ft 


^tt      ,8 


+ 


16fcf 

Sej    ,  12e? 


9e? 


lO/?;6'"1      16/z?6'-1 


9e? 


^f8 
16/??  16/i?   '"        lG/r? 


,ec'-i 


5f'-l  1R^5€<-1 


HVf 


/    ,8 


9*4 


^e 


16hfe'-'  r  16ft3 e,~ 


'     ,8 


r'-i  +  iG/^e'-i  +  iG/.?6 


Finally. 


Part  III. 


r=l        v 


2\2  '  I<t2  =  O,0=/ 


T  T 

1    v-^      o     „g  1    ^~ >  0  ^g 

/=i  -"  /=i 


To  compute  the  moments  of  e,.  let  us  define 


\       _  {e2m       ->(,n-l)  2, 

A(  —  [€t     .  e,  e,  j. 


Following  Engle(19S2).  we  have 


E(e%"  |  $,_! )  =  hm  ]  \(2j  -  1)  =  («0  +  C«?-i )'"  f[(2j  -  U- 


j=i 


3-1 


Expanding  this  expression  establishes  that  the  moment  is  a  linear  combination  of  \t-i 


E{\(  |  $,_2)  =  6 +  .4(6  +  AAt_2). 


Only  powers  of  e  less  than  or  equal  to  2m  are  required,  therefore  A  is  a  upper  triangular 
matrix  and  b  is  a  /??  x  1  vector.  In  general. 


E(X,  |  $,_A.)  =  (/  +  .4  +  .4*  +  ...  +  Ak~[  )b  +  Ak\t-k 
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Because  the  series  starts  indefinitely  far  in  the  past  with  2r  finite  moments,  the  limit  as 
k  goes  to  infinity  exists  if.  and  only  if.  all  the  eigenvalues  of  A  lie  within  the  unit  circle. 
Now, 

lim  E(Xt  |  $,_,)  =  £( A,), 

k  —  oc 

is  an  expression  for  the  stationary  moments  of  the  unconditional  distribution  of  e.  We 
have 

E(\()  =  EE(\,  \$t-1)  =  (I-E(A))-1b. 

Since 

IT      pun-'        /l-3£-(C)-     -Ca0£(C-)V' 
^'EiA))      -{         o  1-£(C2W 

1  —  o  i G  q  o  a  i 

(l-o,)(l-:j(Q'f  +  2a2))         (l-Ql)(l-:Maif  +  2<r2)) 

l-Ol 


() 


(l-Q,)(l-:}(Q?  +  2cr2)) 


and 


/>  =  (3ao,cvo)', 
the  expression  for  the  fourth  and  seconds  moments  are  as  given  in  section  3. 


3G 


